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1.0 INTRODUCTION 

This is the Final Report for Grant NAG 5-325 entitled "Investigation 
of Dynamic Noise Affecting Geodynaraic Instrumentation in a Tethered 
Subsatellite." The author of this report, and the primary analyst, is Dr. 
Gordon Gullahorn. Contributions were also made, and appendices written, by 
Dr. Mario Gross! and Dr. Enrico Lorenzini. The Principal Investigator at 
the inception of the Grant was Dr. Giuseppe Colombo. After Dr. Colombo's 
untimely death. Dr. Gullahorn was made P.I. 

As summarized in the body and appendices of this report, SAO has, 
under this Grant: 


- Formulated a research plan for dynamic noise studies under 
the Grant and a portion of the proposed follow-on. 

- Investigated, in the literature, the current knowledge of 
irregularities in the atmospheric density in the thermosphere. 
Unfortunately, little appears to be known about irregularities 
at the altitudes of interest (120-220 km), the best information 
being available for altitudes around 250-300 km. 

- Calculated effects of the tethered satellite system’s 
internal dynamics on the subsatellite in the frequency range 

of interest. This includes both overall motions (libration and 
attitude oscillations) and internal tether oscillations. 

- Modified the SKYHOOK tether simulation program to 
operate with atmospheric density variations and to output 
quantities of interest. 

- Developed techniques and software for analyzing the results, 
including noise spectral analysis. 

- Begun development of a program for computing a stable 
configuration of a tether system subject to air drag. These 
configurations will be of use as initial conditions for 
SKYHOOK and, through linearized analysis, directly for 
stability and dynamical studies. 

- Performed a simple case study in which the subsatellite traverses 
an atmospheric density enhancement. This has confirmed some 
theoretical calculations, and pointed out some aspects of the 
interaction with the tether system dynamics. 
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Presented, at the first meeting of the Tethered Satellite 
Core Equipment Definition Team at NASA/MSFC in Huntsville, 

Alabama, February 9-10, Some of these results for the purpose 
of defining a core equipment package for the measurement of 
expected acceleration variations in the subsatellite, A copy 
of the presentation by Dr, Grossi is given here as Appendix A. 

Presented, at the NASA Geodynamics Conference (a special session of 
the American Geophysical Union 1984 Spring Meeting) in Cincinnati, 
Ohio, May 14-17, a paper "Feasibility of Gravity Gradient 
Measurements from a Tethered Subsatellite Platform." A copy of 
the presentation by Dr. Gullahorn is included here as Appendix C. 
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2.0 RESEARCH PLAN FOR DYNAMICAL ANALYSIS 

Considering the uncertainties in current knowledge of the fine 
structure of the thermosphere, the effort involved in model development and 
program modification, and the time available for the study, the following 
research plan for the dynamical noise analysis portion of the effort was 
developed. Part of this analysis has been completed under the current 
grant, and part must wait for the proposed follow-on effort. 

Initially, we will consider the subsatellite dynamically as a point 
mass, though it will be subject to appropriate air drag; this is a not 
unreasonable first approximation since the subsatellite is nearly 
spherical. SKYHOOK runs will be made Hith a deployed tether and 
subsatellite moving through a perturbed atmosphere, and the accelerations 
experienced by the subsatellite determined. The subsatellite attitude is 
not explicitly included in SKYHOOK. If deemed appropriate, a simple model 
of the interaction of the subsatellite attitude dynamics and the tether may 
be developed to check for rotational effects. 

It may be of interest to disentangle the direct and indirect effects 
of the atmospheric perturbations. E.g., we might: 

- compute the perturbed accelerations of the subsatellite, 
without actually applying them to the dynamical simulation, 

- apply the perturbed drag only to the subsatellite, but this 
time allow the resultant accelerations of the subsatellite 
to interact with the tether, generating waves, etc., 
which will react back on the subsatellite, 

- apply the perturbed drag only to the tether, to see if the 
indirect effect on the subsatellite is significant. 

These computations should Involve only minor program changes; the major 

effort would be in the analysis of results. This analysis might be useful 
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in developing noise reduction methods. 


The applied atmospheric perturbations will be developed in a stepwise 
fashion; 


- First, a number of simple "units" will be used, suoh as 

- a localized density variation, e.g., a turbule, cell, etc. 

- an extended wave, e.g. a gravity wave 

- a "damped" wave, to simulate those found by Gross and 
Huang (1984; see Section 3.0 below); 

- a step function variation, e.g. a terminator effect. 

Results from these simple models should enhance our understanding of the 
response of the system and the relative importance of different density 
variations. 


- Then, using the best knowledge of the density fluctuations of the 
lower thermosphere, we will generate a sample perturbed atmosphere. 
This will be done by superposing a variety of localized enhancements 
of different scales and strengths. Perhaps several such models may 
be needed to adequately span the possibilities allowed by 

current theory and observation. 

Further modifications that might be desirable in SKYHOOK are: 

- include the subsatellite attitudinal dynamics, including torques 
around the tether attachment point due to air drag; 

- include a model of the gradiometer and its suspension. 

It may also be possible to estimate these effects analytically with less 
effort than modifying and running SKYHOOK. 


Unfortunately, SKYHOOK may not allow fine enough discretization of the 
tether to provide valid results for tether vibrations over the whole 
frequency range of interest. To provide realistic results for periods as 
low as 2 to 3 seconds with SKYHOOK will require so many masses (at least 
twenty) that running simulations will not only be very expensive but have 
inconveniently long turn around. This problem would be aggravated by the 
fact that we will need simulations for long periods of (physical) time in 
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order to get good spectra and to allow the perturbations to fully develop; 
and the fact that the perturbations are constantly varying, even turning on 
and off, which would lead to inefficiency in the SKYHOOK integration 
routine. Besolution of tether effects with 1 to 2 second periods is 
totally infeasible. {This problem is discussed further in Section 8.0.) 

To study system response in this short period regime, we may be able 
to apply the techniques of random structural vibrations (Nigam, 1983) or of 
statistical energy analysis (Lyon, 1975). The structure, in this case, 
would consist of an equilibrium configuration created by the stable 
configuration generator discussed in Section 9.0. This structure is 
subject to random forces, in this case due to the atmospheric density 
fluctuations; the statistics of the induced vibrations are derived from 
those of the forcing function. 
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3.0 ATMOSPHERIC VARIABILITY 

The region in which the subsatellite and tether will be deployed (120 
to 220 km. altitude) is in the lower thermosphere, a region about which 
our direct knowledge appears to be limited. G. Carignan (U. Michigan), J. 
Slowey (SAO), S. Gross (NY Poly) and F. Marcos (AFGL) have generously 
shared their experience and knowledge with us and allowed us to develop a 
brief bibliography. Professor Gross, in particular, sent us copies of 
papers before publication. The relevant papers we have acquired are listed 
in the references. 

There are two forms of possibly relevant data available. The data 
directly on our region of interest appear to come only from perigee 
passages of satellites with elliptical orbits: SPADES (Marcos and 

Champion, 1972), and Atmosphere Explorer-C (Reber, et al. 1975; Potter, 
Kayser and Mauersberger, 1976) . This form of data, brief samples at 
constantly varying altitude, has several disadvantages for our purposes: 
it Is difficult to disentangle the horizontal variations from the vertical; 
thus, the wave vs. random nature of the variations is not convincingly 
demonstrated; and (probably for these reasons) no sophisticated analysis 
has been published. The reason for this form of data is, of course, that 
the atmospheric density is high enough that a satellite would soon decay 
from these altitudes. 

At higher altitudes, still in the thermosphere but substantially above 
the altitudes of interest to us, it is possible to maintain an orbit 
constant for many revolutions. More sophisticated analyses have been made 
of this higher altitude data: for Atmosphere Explorer-E at altitudes near 

300 km (Hoegy, et al., 1979) and for Atmosphere Explorer-C near 250 km 
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(Gross and Huang, 1984; Gross, Reber and Huang, 1984). 

The lower altitude data allow little to be said firmly and 
quantitatively for our purposes (although there are a number of intriguing 
features relating to, for instance, the relative phase of variations in 
different atmospheric components) . Roughly, the variations experienced 
along-trajeetory appear to be in the range of a few percent, up to about 
20%, with distance scales in the range of tens to hundreds of kilometers 
(this is most clearly shown in Figure 3 of Reber, et al.).’ 

The higher altitude data are more complete and better analyzed; 
although their relevance no our purpose may be debated, in some sense it is 
the best we can do. At 250 km altitude. Gross and Huang observe a power 
law spectrum for general variations in the density, and notice some 
discrete waves of finite extent in the horizontal direction: e.g., a wave 

of 250 or 300 km wavelength that extends for only 1250 or 1500 km of a 
total path of 4100 km. They obtain variability over the entire range of 
tens to hundreds of kilometers (about seconds to tens of seconds at a 
velocity of 7.8 km/ sec). Gross, Reber and Huang find variability extending 
to thousands of km with an overall power law spectrum but also some very 
distinct, persistent wave structures, possibly global standing waves. 

One form of information which does not seem to be available is the 
vertic al correlation of density fluctuations. Although the primary 
influence of density fluctuations is certainly the direct drag fluctuation 
on the subsatellite, the effect on the full tether system will depend on 
whether the tether as a whole experiences the fluctuations simultaneously, 
or whether the different regions are perturbed independently. One model 
which could be used to derive the correlation would be to assume that 
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density fluctuations are essentially due to vertical displacement: an 

enhancement, say, seen as the satellite moves horizontally at a given 
altitude would really be due to its seeing air which has been moved up the 
appropriate distance; thus, a higher tether segment would see the nominal 
density from the same distance below its true altitude. This corresponds 
to a gravity wave model, and there is some observational evidence 
corroborating it, at least for wave-like motions. Whether it should apply 
to random (turbulent-like) variations is less clear. 


V. 
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4.0 UNFOBCED SYSTEM DYNAMICS 

It is instructive to look at the dynamics of the system when it is not 
subject to the random drag forcing; we do, however, include the long 
period librations induced by winds in a non-equatorial orbit, as they have 
an Indirect effect on attitude dynamics. Several aspects of the problem 
have been decoupled and analyzed separately; the results are in the 
subsections below. A full self-consistent treatment would be lengthy to 
carry out, though perhaps worth-while; some of the complications are 
discussed in Section 4.2,1. 

4.1 Subsatellite Potation Dynamics; 

The overall librational motion of the tether system is long period 

(orbital or some fraction thereof) and of no direct concern. However, the 

subsatellite also oscillates as a free pendulum with the tether's end as an 

attachment point. Dr. E. Lorenzini has estimated the librational and 

attitudinal oscillations; his calculations are in Appendix B. The period 

of attitudinal oscillation is about 5 seconds, and the associated 

-5 

acceleration about 10 g, (This latter value was computed assuming the 
amplitude of the subsatellite oscillation was the same, in degrees angle, 
as the librational motion.) 

4.2 Tether Elasticity Effects; 

The tether elasticity leads to three sets of vibrational modes In the 
tethered satellite system; 

- spring-mass mode, the subsatellite oscillating as if on 
the end of a massless spring; 
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- longitudinal (stretching) oscillations of the tether; 

- latitudinal ("plucked string") oscillations of 
the tether. 

The latter two are here calculated as if for a tether with fixed ends 
("infinite mass" satellites). The frequencies of the modes are fixed for 
any expected amplitude (so long as the forces are linear), but the 
amplitudes will depend on the specific excitation experienced. The modal 
frequencies are calculated in Section 4.2.1, and a specific excitation 
mechanism is considered in Section 4.2.2. 

4.2.1 Calculation of the Modes: 

Spring/Mass Mode: 

Here we consider the system composed of a large (infinite mass) 

Shuttle, a massless tether with spring constant k = EA/L , and a 

12 

subsatellite of mass m; E - elasticity = 0.7 x 10 for Kevlar, A is the 
tether cross-sectional area, and L the tether length. Then the period of 
oscillation is 


P = 2ji /mL/EA 1 
ms 

which evaluates to = 95 s for our specific mission. 

ms 

Tether Oscillation Modes, Longitudinal and Latitudinal: 

Both modes of oscillation are governed by the wave equation with fixed 


ends: 
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W[“] = u tt - e 2 u n = 0 
u(0,t) = u(L,t) = 0 

where u(x,t) represents the displacement at time t of a physical point on 
the tether whose position at equilibrium would be x. The only difference 
is in the value of c (the speed of propagation) and the interpretation of 
"displacement" (along the tether for longitudinal waves, perpendicular for 
latitudinal waves) . The solution of this equation is a simple exercise in 
separation of variables (see Kaplan, 1981, or almost any applied 
mathematics text) ♦ Any solution will be a linear combination of the 
eigenfunctions: 

u (x,t) - sin(mr x/L) exp (± i w t) 

n 

(s> n = nirc/L 

and these represent the normal modes of the system. Here we are interested 
only in the eigenfrequencies <o n , or the associated periods 


The propagation speeds are given by 

e = /T/fi 1 latitudinal oscillations 
c = /E/p * longitudinal oscillations 




V 
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where T = the equilibrium tension, 
p = the volume density 
p = the linear tether density, = pA. 

12 —2 —3 

For Kevlar (E = 0.7 x 10 dyne cm , p = 1.5 g cm ) and at an altitude 
of 250 km (assuming the equilibrium tension is due solely to linear gravity 
gradient force), these evaluate to 

e = 0.22 km/ sec latitudinal 

c = 6.8 km/ sec longitudinal, 

leading to oscillation periods of 

P « 910/n seconds latitudinal 

n 

P q - 29/n seconds longitudinal. 

Note that both types of waves have modes in the regime of interest, 

from 5 or 10 to 200 seconds. The longitudinal vibrations will have several 

modes near the low period end. The latitudinal oscillations will 

have numerous modes throughout the range, though one may expect that 

the higher modes (short period) will be less strongly excited and more 

highly damped and thus be less distinct. 

The model used in this section is an idealization. It does not 
include the effeots of: 

- damping forces in the tether and elsewhere; 

- finite masses of the end masses (satellite and shuttle); 

- attitudinal dynamics in the end masses; 

- the three dimensional nature of the possible oscillations; 

- the overall (equilibrium) curvature of the tether due to 


_ . ; rf* -V ^ 
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air drag; 

- external forces on the tether, such as gravity gradient or 
air drag (though these will possibly drop out when considering 
small displacements about equilibrium) . 

A full, self-consistent calculation including the above effects should 

perhaps be made. However, it promises to be a lengthy process, including 

coupled sets of ordinary and partial differential equations, and probably 

eventually requiring numerical calculation of the modes. An approximate, 

non-self-consistent calculation of one effect of subsatellite attitude 

dynamics is given in the next section. 


4.2.2 Excitation of Longitudinal Modes by Subsatellite Attitude 
Oscillation: 


To explore how one might get energy into the above modes, we performed 
a sample calculation using the results of sections 4.1 and 4.2.1: 


- Take the subsatellite attitude oscillation period and 
amplitude as given in section 4.1; 

- regard these as a fixed forcing oscillation 

of the tether r s end, and calculate the response in the form 
of longitudinal oscillations; 

- using the derived formulae, compute the tension, and hence 
accelerations, at the subsatellite end of the tether; 
this response oocurs at the forcing frequency and 

at the tether natural mode frequencies. 

We should emphasize that, complicated as the expressions derived may seem, 
this calculation has several limitations: 


- It is not self-consistent; no attempt has been made to 
include the effects of tether oscillations on the sub- 
satellite motion. The small accelerations found give some 
hope that this is not a significant omission except near 
resonances. 

- Only longitudinal waves were considered. The latitudinal ease 
should also be considered; the calculation of the tether response 
would be similar, although calculation of the back-reaction on 
the subsatellite will differ. 

- Resonance effects are not considered. However, the forcing period 
is in the region spanned by the natural frequencies and changes in 




UT ^5^* - 
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mission parameters (or effects not in the calculation) could easily 
put the system close to resonance. Resonance with a high 
latitudinal mode will be even more likely. This clearly 
deserves further study. 


Using the notation of section 4.1.1, the governing equations for the 
tether motion are: 


WEu] = 0 
u(0,t) - 0 

u(L,t) = A sin (to f t); 

a complete problem also requires some form of initial conditions, typically 
specifying u(x,0) and u fc (x,0) at time t=0. 

The calculations are straightforward but lengthy, and we shall only 
give the results. Following Kaplan (1981, pp. 422-423), we try a solution 
of the form 


u(x,t) = Ch I{ .(t)sin(knx/L)] + (x/L) A sintco^ t) 
k=l 

where the h^ (t) are to be determined from the equations. The resulting 
solutions are 


h k (t) = Bsin^t+tf) + (-l) k+1 



r ^ i 

/2A\ 

\k rtj 

*2 ^2 
( 0 . - 6 )^ 
u k £ J 


sin(w f t) 


where B and are constants to be determined (for each mode lc) by the 
initial conditions imposed. To give a concrete example, we set u - u. =» 0 
at t=0. Then, 




ORIGINAL PAGE IS 
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h k (t> - (-1V 


» ft] 


[ (Wjj/Wj, ) sin© k t-3inw f t] 


where © k = knc/L 


= given forcing frequency 

(subsatellite attitude oscillation frequency) 


The tension in the tether is then 


T = Eo u (x,t) 

Jt 


and evaluating this at x = L gives 


[L/CAEo) T = sin(to^t) 


f'-t ft! •£[: 

k=l k f L 


2©, 

©„ 

k 

f j 

tu^ - 

©2 

k 

f j 


We may divide this into the components: 
at the forcing frequency 

T k at the normal mode frequencies © k = knc/L, k = 1,2. 
This tension will have the following effect on the subsatellite: 


- Linear acceleration in the direction along the tangent 
to the tether; broken into the above modes: 

S (> * T () /M 

g 

where M is the subsatellite mass, 0.5 x 10 grams. 


sin©, t 
k 
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From the discussion of Sr. Lorenzini (Appendix B) , He know that the 
forcing frequency corresponds to a period of about 5 seconds. The mean 
angle of the subsatellite to the tether is 1° and the oscillation amplitude 
also about 1°. The amplitude of the forcing oscillations will then be 
about 


A = (75 cm) (1 - cos2°)/2 = .023 cm. 


A brief program was written to evaluate the resultant accelerations. For 
forcing periods near 5 seconds, we note 


- one or two modes near resonance will have 

linear acceleration about 2 to 5 miero-g (10 6 gravities) 

- the rest of the natural modes will have 

linear acceleration less than about 0.5 micro-g, 
decreasing as 1/k for the higher modes 

- the back-reaction at the forcing frequency is roughly 

linear acceleration about 2 to 5 micro-g 


Although these studies have given some definite numbers to use for 
design purposes, the connection to the actual experiment of measuring 
gravity gradients is not clearly defined. This will depend strongly on the 
instrument used and, especially, its suspension system. This connection 
deserves careful study. 
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5.0 SOFTWARE DEVELOPMENT — SKYHOOK 

SAO has developed and maintains SKYHOOK, a general tether simulation 
program. SKYHOOK models the tether as a series of discrete masses 
connected by massless tether segments. The equations of motion are 
integrated in Cartesian coordinates. Most known forces may be taken into 
account, and new versions generated for particular studies. 

A version of SKYHOOK has been prepared which allows superposition of 
density variations on the smooth ambient atmospheric density. Output is 
made to the informational print file, and to a machine readable file for 
plotting and analysis. The non-gravitational acceleration of the 
subsatellite is computed at each output interval, and expressed in two 
orthogonal coordinate systems for output: in one, the axes are defined by 

the local vertical and the orbital motion; in the other, the tangent to 
the tether and the orbital motion are used. The print file also includes 
such quantities as height, ambient unperturbed density and the ratio 
actual/unperturbed for each mass point. This version does not include 
subsatellite attitude dynamics, so rotational velocities are not output. 

The organization has purposely been made general enough that most 
specific density variation models can be accomodated with only minor 
reprogramming, at a limited number of points in the code. New models will 
be added to the old ones, incrementally generalizing the program rather 
than proliferating a sequence of narrowly specific program versions. We 
have developed a set of lexical routines which will allow input parameter 
specification in a "user friendly" flexible fashion, and a set of 
programmer protocols for addition of model enhancements to SKYHOOK. These 
should considerably simplify future changes. 
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The specific perturbation model which has been incorporated so far 
consists of an elliptical region in which the density is enhanced or 
reduced. The perturbation is two-dimensional, in the orbit plane, since 
the distance scales of interest (tens of kilometers or mo-'s) are much 
greater than anticipated out-of-plane motions. As parameters one may 
input: the size of the region (horizontal and vertical axes); the center 

placement (height, and distance "along orbit" from the start of the run); 
and the density ratio perturbed/ ambient. 

A long-standing problem leading to occasional mysterious failures was 
diagnosed and fixed. The atmospheric temperature routine did not always 
handle angle variables correctly; the resultant temperatures could be 
outside the domain of the lookup table used to get atmospheric density, and 
no checks were made to guard against this, so that the contents of some 
random position in memory would become the "density." This led to 
unpredictable results. 

Results of sample simulations are given in section 7 below. 
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6.0 SOFTWARE DEVELOPMENT — ANALYSIS PROGRAMS 

Once a SKYHOOK run has been made and a file of accelerations has been 
produced, we must do something with this data. Several service routines 
have been written to help examine the results. 

6.1 Spectral Analysis 

Gradiometer, accelerometer and other instrument sensitivity are 
commonly quoted as functions of frequency. Thus, the spectral 
decomposition of a sample acceleration profile is, at the least, a useful 
tool for estimating instrument response. 

Experiment with available packaged routines for spectral analysis 
(from IMSL) produced confusing results in test programs, so a routine was 
developed to produce simple spectral results (using, however, the FFT 
routine from IMSL, which performs the clearly defined task of computing a 
discrete Fourier transform (DFT) ) . We describe very briefly the capacities 
of the routines 

If we have N samples of some quantity at times t = 0, dt, 2dt, 3dt, 

... (N-l) dt, then the DFT will produce l+N/2 results at frequencies f = 

0, df, 2df, ... (N/2)df. Here we assume N Is even; a power of 2 (e.g. 

256, 512, 1024, . . .) makes for computational efficiency, but is not 
necessary. The frequencies f are in Hertz. The frequency resolution 
(lowest discernible frequency or longest period) depends on the total 
length of the observed time Interval: df = 1 / (Ndt). The highest 

frequency (shortest period) available depends only on the sampling interval 
dt: (N/2)df = 1 / (2dt) (the Nyquist frequency) . 
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Before we compute the DPT of the sample, we perform a process known as 
"Hanning" to reduce the effects of the finite length of the data stream. 
This consists of multiplying the data by a function that reduces smoothly 
to zero near the edges of the data span ("window"), reducing ringing or 
"leakage" in the Fourier transform domain (see, e.g., Brigham 1974, pp. 
140-146) . Some simple tests we ran confirm that this technique 
dramatically improves the resolution of edges and discrete frequency 
components not at an exact value k.df of the discrete spectral domain. 
Sample data created with an exact frequency k.df (e.g. a simple sine 
function) had their spectra slightly smeared, over the two adjacent 
frequency samples. This seems a small price to pay for the generally 
improved performance, but we should be aware of the decreased resolution; 
in particular, a peak two or four frequency samples wide may actually be 
quite sharp. 

A program was written which reads the file of accelerations produced 
by SKYHOOK. The user may then request the spectrum (amplitude of the 
Fourier transform; this is the square root of the power spectral density) 
for any time interval of data, and the spectra of all six acceleration 
components (three components in each of the two coordinate systems defined 
in section 5.0) are output to a file for printing, plotting or analysis. 
Each component's spectrum is normalized to have a maximum of 1. 

6.2 Other Routines: 

Programs to do printer plots of the accelerations and of the spectra 
from the program in section 6.1 were written. Only the two acceleration 
components tangent to the tether ("along") and orthogonal (in the orbital 
plane; our runs so far have had no out-of-plane component) to this are 
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used. The logarithms of the spectral power are plotted. 

A program was written to difference two of the acceleration files 
output from SKYHOOK. This allowed us to look at only the effects of 
atmospheric perturbations by making identical runs with and without the 
perturbation and subtracting the results. 
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7.0 A CASE STUDY 

A case was run to examine the effects of a single enhanced density 
region. The orbit Has equatorial* with all motion in the orbital plane. 
Ten masses were used: the shuttle* the subsatellite* and eight tether 

masses. Because our standard initial condition generator (D UMBEL) 
considers only the gravity gradient force* we put heavy damping in the 
tether and ran the DUMBEL produced configuration until most of the 
longitudinal oscillations had damped out. The result of this run was used 
as input to the simulations described here* with the damping removed. 

We used a 100 ton shuttle* an 0.5 ton subsatellite deployed downward, 

12 —2 

and a 100 km tether of 0.2 cm diameter Kevlar (E = 0.7 x 10 dyne cm , p 

_3 

- 1.5 g cm ). The Shuttle was in a 220 km altitude, circular, equatorial 
orbit . 


The atmospheric density perturbation consisted of a 20 km radius 
circular region enhanced by 20% over ambient density* with its center at 
120 km altitude (the same height as the subsatellite) and 100 km along 
orbit from the start of the run. 

Two runs were made, each following the tether for 700 (physical) 
seconds. The first run, in which the density perturbation was turned off, 
was for comparison purposes and took 15 minutes of computer time. The 
second run, with the perturbation on, took 19 minutes. 
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7,1 Analysis of the Hesults: 

Since the damped SKYHOOK run which produced our initial conditions had 
not completely relaxed, even the run without the density perturbation 
exhibited substantial oscillations. To isolate the effects of the 
perturbation, we differenced the acceleration output of the two runs and 
worked with the resulting file. In the regime of small disturbances 
explored here, either residual oscillations from "deployment" or due to 
atmospheric perturbations, we expect linearity; thus this scheme should be 
equivalent to starting with a fully relaxed initial state. 

The clearest results were seen looking at the components (1) along the 
tangent to the tether (nearly vertical) , and (2) in the perpendicular 
direction (nearly horizontal, parallel to the drag force) . This choice of 
coordinates most clearly distinguishes the various effects. 

First we consider the accelerations produced by the density 
enhancement In the time domain; we will consider the frequency domain 
below. 

The obvious immediate effect Is an acceleration in the "across tether" 
(horizontal) direction of -0.084 gal (1 gal = 1 cm sec -2 ). This is very 
nearly constant throughout the period when the enhancement is on. The 
pust-enhancement response is limited to the range [-8.1 micro-gal , +49 
micro-gal), i.e. down by more than three orders of magnitude, A printer 
plot of this post-enhancement response is given in Appendix D, Figure D-l, 
The overall variation Is seen to be quite smooth, with little short scale 
oscillation. The Immediate response during the enhancement is not plotted, 
since the post-enhancement response would not then be discernable. 
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The entire acceleration response in the "along-tether" direction 
(vertical) is shown in Figure D-2. The total range of response is 1-0.086 
gal, +0.096 gal]. During the enhancement, the acceleration along-tether 
rises gradually to 0.024 gal. This brief perturbation has excited the 
tether, and in bhe post-enhancement response we see substantially larger 
accelerations produced by the back-reaction of the tether rather than 
directly by the enhancement. Note that these post-enhancement 
accelerations in the along-tether direction are much larger than the 
across-tether accelerations. Also, the along-tether response shows much 
more variation on short time scales, a fact which will also be apparent 
from the spectral analysis below. 

The spectra of both components were computed as outlined in Section 
6.1. Data for 512 points (at one second intervals) immediately 
post-enhancement were used, giving a frequency resolution of about .002 Hz. 
This is actually somewhat misleading since the Hanning window used to 
suppress "leakage" from sharp spectral peaks also smears them by about 1 
frequency resolution element on either side; thus the true resolution is 
closer to .006 Hz. Printer plots for the across-tether and along-tether 
components are given in Appendix D, Figures D-3 and D-4 respectively. 

These are the logarithms of the spectral amplitude. A portion of the 
along-tether spectral decomposition, labeled so as to be the logarithm of 
the power spectral density, is shown in Figure 7-1. 


LOGm (POWER) 




ACCELERATION IMPULSE RESPONSE P.S.D, 


LONGITUDINAL MODES 




FREQUENCY (HZ) 


Figure 7-1. Power Spectral Densidy of the subsatellite acceleration after encountering 

an atmospheric density enhancement (see discussion in text). The along-tether 
component is plotted; for a more extensive plot and for the perpendicular 
component, see Appendix D. The mass-spring mode ("MS") and the expected 
eight longitudinal oscillation modes are marked. 


original 13 
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Comparing Figures D-3 and D-4, we see directly that the aeross-tether 
component shows only a few indistinct peaks at low frequency (long period) , 
while the along-tether component shows several strong peaks down to a 
period of about 5 seconds, before decreasing to a random tail. The 
physical nature of the peaks is labeled in Figure 7-1. One is seen to be 
the spring-mass mode discussed in Section 4.2.1, with a period of about 100 
seconds, as close to the predicted value of 95 seconds as can be 
determined. We also see eight modes of longitudinal tether oscillation; 
their frequencies, and the frequencies predicted from the model of the 


continuous tether with fixed 

ends, are: 

Mode 

Theoretical 

Simulation 

1 

0.034 Hz 

0.037 Hz 

2 

0.068 

0.068 

3 

0.102 

0.100 

4 

0.137 

0.127 

5 

0.171 

0.148 

6 

0.205 

0.170 

7 

0.239 

0.184 

8 

0.273 

0.193 


We see that the agreement is quite good for the lower modes, and gradually 
deteriorates for the higher modes. The disagreement i 3 probably due to the 
discretization (point masses and springs model) used in the simulation, 
although the limits of the theoretical model (primarily the assumption of 
fixed ends) also introduces some error in the theoretical figures. Note 
that precisely eight modes are expected for the simulation, as observed: 
ten masses were used, eight internal plus the shuttle and subsatellite, and 
one will expect one longitudinal mode for each internal mass. The 
latitudinal modes do not show up at all in the simulation results. This is 
probably due to their low frequency: the low order, stronger, modes will 

be confused with the spring-mass mode and all the modes will be so closely 
spaced as to be unresolved from our data. 
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8.0 ADEQUACY OF SKYHOOK FOR VIBRATION STUDIES 

The above case study, with ten masses run for twelve minutes of 
simulated time, and with only one density enhancement, was a significant 
computational effort. But as we can see from the table given In Section 
7.1, the highest frequency adequately resolved is 0.15 to 0.20 Hz. 

More complete results will require accurate simulation of frequencies 
to 0.5 or even 1.0 Hz. We can compute two estimates of the number of 
masses needed: the first is based on the naive assumption that all the 

theoretical modes are the same as the simulated modes; and the second on 
the observation that about half the simulated modes adequately approximate 
the real (theoretical) modes. These estimates give 17 masses and 32 masses 
respectively to achieve 0.5 Hz resolution; and 32 and 61 masses for 1.0 
Hz. 


Additionally, we will wish to make runs for much longer periods to 
allow the full effects of the noise to build up, and to include the 
dynamics of the subsatellite and the gravimeter suspension. A realistic 
atmospheric model will include many more perturbations, each of whioh will 
cause the integrator used in SKYHOOK to have at least some difficulty in 
continuing and will increase the time required for solution. 

Experience has shown that SKYHOOK simulations for more than ten masses 
become expensive very rapidly, as well as leading to inconvenient turn 
around time. Twenty masses is about the limit feasible even, in a very 
smoothly evolving case, and thirty to sixty mass runs would require a large 
dedicated machine unless the programs efficiency can be very substantially 
enhanced. 




, aSk JSi * 
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One non-SKYHOOK option which is being explored is the use of 
statistical methods for dynamical systems, such as random vibrations 
(Nigam.1983) or statistical energy analysis (Lyon, 1975). In broad 
outline, these methods require a model of the system in an equilibrium 
state, and a statistical description of the perturbing forces (e.g. 
atmospheric density irregularities) . As output, one achieves a statistical 
description of the system response. This is used to derive quantities of 
engineering interest. Some steps have already been taken, for another 
purpose, which would help implementation (Section 9.0 below), though 
significant study would be required before useful results could be 


obtained 
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9.0 SOFTWARE DEVELOPMENT -- EQUILIBRIUM CONFIGURATION GENERATOR 

We have partially developed a program to generate equilibrium deployed 
configurations with large numbers of masses and taking into account the 
non-perturbed atmospheric drag. Other forces could be easily added if 
deemed advisable. 

The original motivation was to generate initial conditions for use 
with SKYHOOK. The current program used to generate initial configurations, 
DUMBEL, does not account for drag forces and has other limitations. Thus, 
the first stages of the SKYHOOK run consist of a long period spent damping 
vibrations before any significant analysis of the perturbed response can be 
made. This i3 expensive enough in both computer and person time, that an 
effort to generate usable initial configurations already in equilibrium is 
justified. 

An additional benefit is that the equilibrium configurations thus 
generated could be used as the systems for the statistical perturbation 
analyses mentioned in Section 9.0. In addition to the configuration we 
would need to compute the linearized force matrix, but this would be a 
simple addition to the program. 

The linearized system could also be used in stability analyses by 
simply computing the eigenvalues of the force matrix using available 
subroutine packages. 

The program has been fully designed in such a manner that convergence 
to a true solution seems assured from both a physical and an algorithmic 
viewpoint. Schematically, the program consists of several levels of 
subroutines each of which solves for a given variable or set of variables 
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by an iterative or exact process; when iterative processes are used, the 
lower level subroutine may be called upon many times to solve its own 
sub-problem with slightly different pi 1 vneters. Some of the inner loops 
have been coded (a conscious departure from standard top-down programming 
practice), and completion should take about one week. 


V 
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10.0 SUGGESTIONS FOR SYSTEM DEVELOPMENT 

From the above analyses and simulations we see that the acceleration 
perturbations experienced by instruments in an atmospheric mission 
subsatellite are certain to be substantial and complex. Elimination or 
compensation for these perturbations will be crucial to the success of a 
gradiometer mission. 

It is anticipated that the longitudinal tether vibrations excited by 
the atmospheric density variations can be damped either by physical devices 
(appropriately chosen spring-dashpot combinations) or by manipulation of 
the tether reel on the deployer as discussed in Colombo, et al A (1984, 
section 5.1) in a different context. 

It is likely that measurements of the atmospheric density 
perturbations will allow removal of some effects in post-processing. This 
measurement could be made either with an external instrument or simply by 
observing the varying drag with an on board accelerometer. 

Rotations, as well as linear acceleration, will be induced in the 
subsatellite. These effects should be more clearly understood. Rotation 
of the gradiometer instrument will give false readings. However, If all 
nine gravity gradient tensor components are measured, it is In theory 
possible to measure the rotation and remove its effects leaving only the 
actual gravity gradient. 

In any event, the suspension of the gradiometer within the 
subsatellite is crucial. It cannot be simply attached rigidly to the 
instrument platform. The gradiometer itself will exercise a significant 
"common-mode" rejection of linear acceleration noise: that is, since all 





accelerometer components of the gradiometer are perturbed by the 3ame 
amount, the measured gradient (the difference of accelerations) is not, in 
theory, perturbed. The degree of this rejection achieved by various 
instruments must be determined. More difficult to circumvent are the 
effects of rotation. Even if the subsatellite itself does not rotate, a 
one-point suspension would cause gradiometer rotation in response to linear 
accelerations. A suspension must be designed to obviate this problem. 
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A Presentation Made By 
Dr. Mario D. Grossi 
to 

Tethered Satellite Core Equipment Definition Team 
NASA/MSFC Huntsville, Alabama 
February 9«10, 1984 


Core Equipment for TSS: 
Sensor of Dynamic Noise 
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CORE EQUIPMENT FOR T.S.S, 


SENSOR OF DYNAMIC NOISE 


FEBRUARY 1984 

SMITHSONIAN INSTITUTION 
ASTROPHYSICAL OBSERVATORY 
CAMBRIDGE; MASSACHUSETTS 02138 


RECOMMENDED SENSOR SYSTEM 


• THREE LINEAR ACCELEROMETERS, ONE FOR EACH OF THE 
THREE X, Y, Z AXES 


• THREE GYROS, EACH WITH AXIS COINCIDENT WITH ONE 
OF THE THREE X, Y, Z AXES 


• A SET OF TRANSLATIONAL AND ROTATIONAL DAMPERS 10 
ISOLATE TETHERED SUBSATELLITE FROM ACCELERATION 
NOISE CAUSED BY TETHER DYNAMICS 




TYPICAL USERS' DEMANDS 


EXAMPLE OF USERS : GEODYNAMICS COMMUNITY, WITH GRAVITY GRADIOMETER AS INSTRUMENT ON-BOARD TETHERED 

SUBSATELLITES 


(A) REQUIREMENTS FOR LIMITATION IN UNCORRELATED NOISE, PROOF MASS - TO - PROOF MASS , 
VALUES COMPUTED FOR 1 SEC INTEGRATION TIME, GRADIOMETER BASELINE 50 CM *: 


0.35 

-9 -2 

10 cm sec = 0.35 

-1° „ 

10- 12 

-13 

g 

to achieve a threshold 

sensitivity of 

10" 2 

10" 3 

-4 

EtJtvBs 

Units 

0.35 

10 " = 0.35 

-11 

10 

-14 

g 

11 

II 

11 

11 
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0.35 

10 " = 0.35 

10 

g 

II 

II 

* 11 

10 
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(LINEAR ACCELERATION NOISE) 


BECAUSE AT THE VERY BEST WE WILL HAVE A THRESHOLD SENSITIVITY OF 10 g / ( MAYBE WE 

CAN HAVE A THRESHLOD SENSITIVITY AS GOOD AS 10 -1 0 g /pzJ,'WE HAVE TO RELY IN THE COMMON-MODE 
REJECTION MECHANISM ( MASS-TO-MASS NOISE CORRELATION) TO' ACHIEVE THE FOLLOWING DEGREES OF 
NOISE CANCELLATION : 

IF USING A 10 _8 g /l/H? IF USING A 10~ 10 g /||hz* GRADIOMETER 

ACCELEROMETER ACCELEROMETER SENSITIVITY GOAL 


A FACTOR OF 

3 

10 4 

A FACTOR OF 

3 

CM 

O 

i—l 

io " 2 

EU 

ti 

3 

10 5 

t! 

3 

10 3 

10~ 3 

ft 

T? 

3 

10 6 

tt 

3 

10 4 

10" 4 

Tt 


THESE SPECIFICATIONS ARE DIFFICULT TO MEET WITH A SIMPLE AND LOW-COST ACCELEROMETER. 


TYPICAL USERS' DEMANDS (cont.) 


WE COULD CONSIDER USING SUCH HIGH-PERFORMANCE ACCELEROMETERS AS THE ONES BEING 
DEVELOPED IN FRANCE BY ONERA, MATRAS, ETC. [THRESHOLD SENSITIVITY 10-13 


OUR RECOMMENDATIONS ARE, HOWEVER, THAT WE SHOULD KEEP WITHIN PRACTICAL BOUNDS 
ACCELEROMETER SENSITIVITY. WE SHOULD PERFORM TRADE-OFF ANALYSIS BETWEEN THIS 
SENSITIVITY AND ACHIEVABLE DEGREE OF COMMON-MODE REJECTION, WITH THE AIM OF 
ESTABLISHING A FAIR SHARE OF THE BURDEN (BASED ON FEASIBILITY, PRACTICALITY 
AND COST) BETWEEN ACCELEROMETER REQUIRED PERFORMANCE AND COMMON-MODE REJECTION 
SPECIFICATIONS. 


AT THIS POINT, WE EXPECT THAT 10“8 g/Sz (POSSIBLY 10~10 g/i/ftl) SHOULD BE THE 
SENSITIVITY TO BE CONSIDERED IN THE PRE- PROCUREMENT SEARCH FOR A SUITABLE 
ACCELEROMETER. 


ANALYSIS OF PRACTICALLY ACHIEVABLE COMMON-MODE REJECTION IS, AT PRESENT, AN 
URGENT, OPEN REQUIREMENT. 

AMOUNT OF ACHIEVABLE REJECTION: TBD 

(B) REQUIREMENTS FOR LIMITATIONS IN ROTATIONAL ACCELERATION NOISE, VALUES OF GYRO 
ANGULAR VELOCITY STABILITY REQUIRED IN 1 SEC INTEGRATION TIME: 

Acu < 3.16 10 -6 rad/sec 

Aid < 10 - 6 rad/sec 

Ato < 3.16 10"? rad/sec 


to achieve a threshold sensitivity of 10 - 2 Eotvos Units 
to achieve a threshold sensitivity of 10-3 Eotvos Units 
to achieve a threshold sensitivity of 10“4 Eotvos Units 


A- 5 


3 


TYPICAL USERS' DEMANDS (cont.) 


WITH A GYRO THAT HAS A VELOCITY STABILITY OF I MILLIDEGREE/HOUR, WE HAVE 


A 10 1 _ Q 

tO = #. = 5 10 5 rad/sec 

57 3600 

THIS AMPLY SATISFIES REQUIREMENTS SET FORTH IN PREVIOUS SLIDE. 

WITH A GYRO THAT HAS A VELOCITY STABILITY OF 3 MILLIDEGEEE/HOUR WE HAVE 
Ad>£' 10-8 ra d/sec. THIS GYRO ALSO SATISFIES ( WITH A SMALLER MARGIN, 
THOUGH) THE REQUIREMENTS ABOVE. 


WE CONCLUDE THAT STATE-OF-THE-ART GYROS HAVE A PERFORMANCE GOOD ENOUGH 
TO BE USED IN T.S.S. DYNAMIC NOISE SENSING INSTRUMENTATION. NO PROBLEMS 
HERE. 


(C) SPECTRAL SENSITIVITY OF INTEREST TO GEODYNAMICS : THE BAND BETWEEN 

10 HZ AND 1 HZ ( IF IT WOULD MAKE A DIFFERENCE FOR INSTRUMENT 
COMPLEXITY, THIS BAND COULD BE REDUCED TO 10 " 2 HZ-to-0.2 HZ). 
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LIST OF SOURCES OF DYNAMIC HOISE 


o ATMOSPHERIC DENSITY GRANULARITIES, AS A CAUSE OF TIME -VARIABLE DRAG 

o TETHER'S PENDULAR LIBKATIONS, IN-PLANE AND OUT-OF-PLANE 

30 NATURAL MODES OF LONGITUDINAL (STRETCHING) OSCILLATIONS (TENSION WAVE MODES) 

o ACCELERATIONS IN PAYLOAD MODULE OF SUBSATELLITE DUE TO ITS ROTATIONAL DYNAMICS 
WITH RESPECT TO TETHER 

o LONGITUDINAL OSCILLATIONS DDE TO TETHER CONTROL 
o LATITUDINAL OSCILLATIONS 

o GROSS MOTION OF SUBSATELLITE ALONG DIRECTION OF TETHER ( SPRING- MASS MODEL) 
o FLUCTUATIONS OF ARTIFICIAL GRAVITY ON-BOARD SUBSATELLITE 


CHARACTERIZATION OF EXPECTED DYNAMIC NOISE 


o ATMOSPHERIC DENSITY GRANULARITIES 

( WE ASSUME 10 % ARIR DENSITY SPATIAL FLUCTUATIONS) 


100 Kin TETHER, DOWNWARDS, 
SHUTTLE AT 230 KM, 28° 


INCLINATION (CASE A) 


20 Km TETHER, UPWARDS, 
iUTTLE AT 295 KM,. 28° 
NCLINATION (CASE B) 


o TETHER PENDULAR LIBRATIONS 


CASE A 


ACCELERATION 

MAGNITUDE 


FREQUENCY 

SPECTRUM 


If' 



0.5 10~ 4 g 

4 10“ 3 HZ to 0.25 HZ 


-7 

, 


0.5 10 g 

same as above 


-3 

. *' » -4 


10 g (IN PLANE) 

freq. « 4 10 HZ 


-5 

-4 


10 g (OUT OF PLANE) 

" 2 10 HZ 

l 

max 

_ c 

=4 10 rad/sec 

freq. = 4 10 HZ 

1 

max 

-5 

=■2 10 rad/sec 

" 2 10~ 4 HZ 






— * 


2 , 


CHARACTERIZATION OF EXPECTED DYNAMIC NOISE (cont.) 


o TETHER PENDULAR LIBRATIOHS (CONT. ) 


CASE B 


ACCELERATION 

MAGNITUDE 


FREQUENCY 

SPECTRUM 


3 

2 



-5 

g ( IN PLANE) 


-4 


10 

freq. 

= 3 10 

HZ 

, -5 



-4 


10 

g (OUT OF PLANE) 

ir 

4 10 

HZ 

« 7 

-6 , 

10 rad /sec 

freq. 

» 3 10 -4 

HZ 

= 4 

10 ^ rad/sec 

U 

4 10~ 4 

HZ 


NOTE : THESE ACCELERATIONS ARE EXPECTED TO BE EFFECTIVELY CANCELLED BY FILTERING 

AND BY THE USE OF TRANSLATIONAL AND ROTATIONAL DAMPERS ON BOARD SUBSATELLITE. 
ESTIMATE OF OVERALL ACHIEVABLE DEGREE OF CANCELLATION : TBD. 


o NATURAL MODES OF LONGITUDINAL (STRETCHING) OSCILLATIONS (TENSION WAVE MODES) 


CASE A 


5.3 10“ 6 g 


O.I66 HZ to 0.73 HZ 


CASE B 


5 10 



HZ to 1.25 HZ 




CHARACTERIZATION OF EXPECTED DYNAMIC NOISE (cont.) 


o ACCELERATIONS XN PAYLOAD MODULE OF SUBSATELLITE DUE TO ROTATIONAL DYNAMICS 
WITH RESPECT TO TETHER 


CASE A 


ACCELERATION 
MAGNITUDE 
-3 


CASE B 


FREQUENCY 

SPECTRUM 





10 

g (IN PLANE) 

freq. » 

0.2 

HZ 




-5 







4 

10 

g (OUT OF PLANE) 

TT 

H 

it 




-5 
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2 

10 

g (CENTRIFUGAL) 

freq. = 

TED 


K, 

inax 

= 

00022 rad /sec 

freq, = 

= 0.2 

HZ 

1 CO 1 


=2 

10 ~ 

1 

rad /sec 
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1 ] {max 
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10 
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10 
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HZ 
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NOTE: THESE ACCELERATIONS ARK EXPECTED TO BE EFFECTIVELY CANCELLED BY FILTERING AND 

BY THE USE OF TRANSLATIONAL AND ROTATIONAL DAMPERS ON BOARD SUBSATELLITE . 
ESTIMATE OF OVERALL ACHIEVABLE DEGREE OF CANCELLATION : TBD . 


CHARACTERIZATION OF EXPECTED DYNAMIC NOISE (cont.) 


4 . 


o LONGITUDINAL OSCILLATIONS DUE TO TETHER CONTROL 


ACCELERATION FREQUENCY 

MAGNITUDE SPECTRUM 


CASE A & CASE B 


TBD 


FREQ. EXPECTED TO BE 
OF THE ORDER OF 10~ 4 HZ 


o LATITUDINAL OSCILLATIONS 


CASE A 


TBD 


freq. =5.4 10 HZ 


CASE B 


TBD 


freq„= 9.4 10 HZ 


o GROSS MOTION OF SUBSATELLITE ALONG DIRECTION OF TETHER ( SPRING- MSS MODEL) 


CASE A 
CASE B 


-3 

2 10 g 

TBD 


freq. = 10“ z HZ 
freq.= 3 10~ 2 HZ 


o FLUCTUATIONS OF ARTIFICIAL GRAVITY ON BOARD SUBSATELLITE 


CASE A 
CASE B 


6.8 10" 5 g 

1.35 10~ 5 g 


freq.= 3.7 10 HZ 

freq. =3.7 10~ 4 HZ 


INITIAL RECOMMENDATIONS FOR 
BASIC PERFORMANCE SPECIFICATIONS 
OF SENSOR SYSTEM 


WE HAVE SEEN IN THE FIRST PART OF THIS PRESENTATION HOW TO FORMULATE SENSORS' SPECIFICATIONS 
FROM THE VIEWPOINT OF TYPICAL POTENTIAL USERS, SUCH AS THE GEODYNAMICS COMMUNITY. WE SAW 
THAT THE SPECIFICATIONS TO BE IMPOSED TO THE GYROS ARE QUITE MILD, WHILE THE ONES FOR THE 
LINEAR ACCELEROMETERS ARE VERY THOUGH. 


IN THE SECOND PART OF THIS PRESENTATION, WE HAVE SEEN ESTIMATES OF THE DYNAMICAL NOISE TO 
BE EXPECTED IN THE T.S.S. SUBSATELLITE FROM TETHER- RELATED SOURCES AS WELL AS FROM 
GEOPHYSICAL SOURCES. 


BY ASSUMING, INITIALLY, THAT WE DO NOT HAVE IN USE MEANS TO ABATE THESE SOURCES (A 
PESSIMISTIC ASSUMPTION) , HERE UNDER ARE OUR RECOMMENDATIONS FOR THE SENSOR SYSTEM 
REQUIRED PERFORMANCE: 


o THRESHOLD SENSITIVITY OF ACCELEROMETERS: nr 8 -to-l£T 10 g/4Iz, IN THE 

FREQUENCY BAND 10-4 Hz-to-1.0 Hz, UPPER BOUND OF DYNAMIC RANGE 5 10-2 
g.ALSO SUITABLE FOR LAB TESTS IN 1 g ENVIRONMENT. 

o THRESHOLD SENSITIVITY OF GYROS (EXPRESSED IN TERMS OF’ REQUIRED ROTATIONAL 
VELOCITY STABILITY) : 1 to 3 MILLIDEGREES/HOUR, IN FREQUENCY BAND 10-4 Kz 

to 1.0 Hz. DYNAMIC RANGE: TBD. 


WE ALSO RECOMMEND THAT TRANSLATIONAL AND ROTATIONAL DAMPING PROVISIONS BE ADDED TO T.S.S. 
SUBSATELLITE TO ISOLATE IT FROM TETHER SYSTEM DYNAMICS. SPECIFICATIONS, AND REVERBERATION 
OF EXTENT OF FEASIBLE ISOLATION IN SPECIFICATIONS ABOVE,: TBD. 


"MINIMUM PERFORMANCE" 
SENSOR SYSTEM 


BY REVIEWING VALUES OF ACCELERATION MAGNITUDES FROM THE VARIOUS SOURCES OF DYNAMIC 
NOISE TAKEN INTO ACCOUNT IN THE PREVIOUS SLIDES, THE FOLLOWING "MINIMUM PERFORMANCE" 
SENSOR SYSTEM CAN BE SPECIFIED, UNDER THE ASSUMPTION THAT FOR THE FIRST T.S.S MISSION 
IT MIGHT BE SUFFICIENT TO MEASURE THE RAW ACCELERATIONS, WITHOUT USING DAMPERS NOR 
FILTERS : 


o 10" 6 g / VHz ACCELEROMETERS IN THE FREQUENCY RAND 10~ 4 to 1.0 HZ, 
WITH DYNAMIC RANGE UPPER BOUND : 5 10 g . Also suitable for lab 

tests in 1 g environment; 

o 3 MI LLI DEGREES/HOUR GYROS, FREQUENCY BAND 10~ 4 TO 1 HZ, DYNAMIC 
RANGE : TBD. 


WITH THIS INSTRUMENTATION , COLLECTED DATA WILL BE OF LIMITED USE TO USER COMMUNITIES. 
HOWEVER, THEY WOULD BE USEFUL IN DESIGNING FOLLOW-ON VERSIONS OF THE INSTRUMENTATION 
PACKAGE WITH PERFORMANCE APPROACHING UPPER BOUNDS STATED IN .PREVIOUS SLIDE. 


OPEN ISSUES 


(CANDIDATE TOPICS FOR SPECIAL STUDIES) 


o Quantitative estimate of noise cancellation to be expected from common-mode 
rejection in typical users' instrument such as a gravity gradiometer 


o Formulate more accurate model of acceleration noise due to atmospheric density 
fl uctuations 

o Definition and design of translational ar.d rotational dampers 
o Quantitative estimate of degree of isolation achievable with dampers above 

o Quantitative estimate of overall degree of isolation by dampers and filters 


o Re-statement of accelerometers' threshold sensitivity requirements, taking into 
account overall achievable degree of noise abatament, inclusive of common-mode 
rejection ? 


Appendix B 


Some Evaluations Related to Dynamic Noise 
in the Payload Module of the Satellite 

by 

Dr. Enrico Lorenzini 


This appendix derives expressions for natural frequencies 
related to the librational and attitudinal dynamics of the tether/ 
subsatellite system. The expressions are evaluated for two cases, 
an atmospheric (downward) probe and an electrodynamic (upward) 
mission. See Figure B.l for a picture of the system and definitions 
of the symbols. 







The attitude free oscillation frequencies of the satellite in the 
orbital plane (0) and in the transverse plane (£) are respectively 
given by: 
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“0 - J 


ifli = 


0 

Tr 


3iA W r 


3£l 2 m 

sat 


(1.1) 
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where T is the tether tension and r the satellite radius. 

For small oscillation angles we can express the three components 
of the accelerations as: 


0sat 


- -b w 2 A 0 sinCto A t + y fi ) in-plane acceleration 


V 


r 0 


a. ^ = -b a: 2 A.sinfdJ.t + y,1 out-of-plane acceleration 
<psat (p (j) <p qr 

a zsat = b[a, 0 A e cosCco e t + y 0 ) + + v^] 2 centrifugal 

3.CC61 S3T3. L 10T1 

Cl. 2) 

where A^ and A^ are the amplitudes of the rotational accelerations. 

The maximum values of the accelerations are given by: 
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a < b [w„A- + to, A.]' 

zsat 1 0 0 d> o 
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With formulae (1.3) we are able to give an estimation of the maximum 


acceleration provided by the satellite rotational dynamics while the 
angular velocities are given by (maximum values) : 
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(1.4) 

•^sat 

max 

“A 

yaw oscillatior. 



The accelerations depend strongly on the mission scenario and on 
the damping of the rotational oscillations of the satellite. Presently 
no damping devices for rotational oscillations are foreseen on the 
satellite so that we will give an estimation of the oscillations in 
this case. With damping devices, on the contrary, these oscillations 
will be damped out. 


Accelerations due to the Librational Dyn a mics of the Tether 

The tethered system, because of the atmosphere or electrodynamic 
drag acting on it, librates both in-plane and out-of-plane with very 
low frequencies. 

Expressions of the tether librations are given by: 


0 = P 0 sinCO Q t + y 0 ) 
♦ = sinitlf + Y(f) ) 


( 2 . 1 ) 


where g 0 and are the amplitudes and £2 Q , £2^ the angular frequencies. 
Maximum angular rates are given by: 


0 = PA 

'max 0 0 


Id) I = P , 
|T| max $ <j> 


( 2 . 2 ) 


Accelerations can be easily derived in the small oscillation case as: 
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Preliminary Evaluation of the Effect of Rotational Dynamics and Tether 
Librations for some Mission Cases 

Let us consider some typical mission cases that give an idea of 
what the acceleration values could be in the case that no rotational 
damping devices are operating on the satellite. 

Satellite Characteristics 

m sat = 500 k * 

J a = J. 100 kg m 2 

6 <p 

r = satellite radius = 0.75 m 

b = distance between satellite c.g. and gradiometer = 0.35 m 

1st Case : Atmospheric Mission (Downward) 

Orbit inclination = 28° 

Tether length - 100 km 

Orbiter altitude = 230 km 

Orbiter rate = SI = 1.178 x 10” 3 rad/sec 
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and their maximum values by: 
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The satellite oscillation amplitude and angular frequencies are 
approximately given by: 


1 ° 

0.04° 


0 , 


= 20 = 2.36 x 10" 3 rad/sec 


Tether libration 
amplitudes and 
forced frequencies 
(atmospheric drag) 


0, * 0 = 1.178 x 10" 3 rad/sec 




A„ “ l c 


0. 04 


Attitude oscillation 
amplitudes and 
free frequencies 




= 1.25 rad/sec (Period = 5 sec) 


o> 


<t> 


= 1.25 rad/sec (Period a 5 sec) 


An estimation of the maximum values of the related accelerations and 
frequencies is therefore as follows: 


| a 0 sat| = °* 010 m/sec 2 ; f Q =0.2 hertz (t = 5 sec) 

1 ’max ' L 


4>sat 


= 3.8 x 10“ 4 m/sec 2 ; f =0.2 hertz (x = 5 sec) 


max 


'4»sat 


a zsat = x 10 ~ 4 m / sec2 5 f 2S at ~ difficult to be estimated 

mav 


'max 


a Q = 0.009 m/sec 2 ; f Q = 3.74 x 10 4 hertz (x = 45 min) 

max 


a, = 9.7 x 10" 5 m/sec 2 ; f = 1.87 x 10“ 4 hertz (x = 90 min) 

^’max ^ 
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and the velocities: 


0 = 0.022 rad/sec • same frequencies as above 

s r max 


r sat 


= 8.73 x lC^rad/sec ; 


max 


|e| = 4.12 x 10“ 5 rad/sec ; 

JTlaX 


|4»| max = 2.06 x 10 5 rad/sec ; 


2nd Case : Electrodynamic Mission (Upward) 


Orbit inclination = 28° 

Tether length = 20 km 

Orbiter Altitude = 295 km 

Orbital rate = 1.175 x 10“ 3 rad/sec 

Current intensity ~ 0.35 Amp 


Satellite oscillation amplitudes and angular frequencies are 


approximately as follows: 

B 0 = 0.2° 

% " °- 10 

= 2.04 x 10“ 3 rad/sec 

t? 

fi, = 2.35 x 10“ 3 rad/sec 

<P 

A 0 * °* 2 ° 

A. = 0.1° 

<P 

ajg = 0.56 rad/sec (Period = 11 sec) 


Amplitudes and 
frequencies of the 
tether librations 


Amplitudes and free 
oscillation frequencies 
of the rotational 
satellite dynamics 


id 




0.S6 rad/sec (Period = 11 sec) 


v 
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Acceleration and frequencies in the payload modules are given by: 
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'max 
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6s at 

= 8.91 x 10” 2 hertz (x = 11 sec) 

cj> sat 

'max 

1.92 x 10“ 4 m/sec 2 ; 

f 

<|>sat 

= 8.91 x 10 2 hertz (x = 11 sec) 

'zsat 

'max 

3.01 x 10" 6 m/sec 2 ; 

f 

zsat 

= difficult to be evaluated 

el 

'max 

2.88 x 10~ 4 m/sec 2 ; 

f e = 

3.25 x 10 -4 hertz (x = 51 min) 

<J) (max 

1.93 x 10” 4 m/sec 2 ; 

f <(> = 

3.74 x 10~ 4 hertz (x = 45 min) 


and the angular velocities: 

|8 


sat 


sat 


max 


max 


'max 


= 0.002 rad/sec j 
= 0.001 rad/sec ; 
- 7.12 x 10~ 6 rad/sec ; 


I ^ I max = 4,1 x 10 " 6 rad/sec > 


same frequenices as above 


The yaw oscillation velocity has not been calculated because it depends 
on many parameters that still require a complete definition. 
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APPENDIX C 


A Presentation Made By 
Dr. Gordon E. Gullahorn 
to 


NASA Geodynamics Conference 
American Geophysical Union 1984 Spring Meeting 
Cincinnati, Ohio 
May 14-17, 1984 


Feasibility of Gravity Gradient Measurements 
from a Tethered Subsatellite 


FEASIBILITY OF GRAVITY GRADIENT MEASUREMENTS 


FROM A 


TETHERED SUBSATELLITE PLATFORM 


G. E, GULLAHORN 
F. FULIGNI* 

M. D. GROSSI 


MAY 1984 


SMITHSONIAN INSTITUTION 
ASTROPHYSICAL OBSERVATORY 
CAMBRIDGE, MASSACHUSETTS 02138 


* VISTING SCIENTIST, IFSI-CNR, FRASCATI, ITALY 
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"TRADITIONAL" GRAVITY MEASUREMENT 


EARTH SURFACE : 

- MEASURE ABSOLUTE ACCELERATION DIRECTLY 

— USE MASS ON SPRING, PENDULUM, FALLING BODY 
— SHORT DISTANCE SCALE EFFECTS 

- HISTORICAL BASE 

SATELLITE ORBIT OBSERVATION : 

- LASER RANGING, DOPPLER RANGING, OPTICAL POSITIONS 
— LONG SCALE EFFECTS 

- INTEGRATED EFFECT OVER MANY ORBITS 

- STATION COORDINATES 

~ POSSIBILITY OF DIRECTLY OBSERVING CONTINENTAL DRIFT 


GRAVITY GRAD I OMETRY 


- ORBITING PLATFORM 

- MEASURE SPACE DERIVATIVES OF GRAVITATIONAL ACCELERATION 

- DIRECT MEASUREMENT 

- MEASURE ANOMALIES AT ALL SCALES DOWN TO ORBIT ALTITUDE 

- ORBIT CAUSES US TO SCAN ACROSS FEATURES 

PROPOSED PLATFORMS : 

- HIGH ALTITUDE FREE SATELLITE 

- LOW ALTITUDE "DRAG-FREE" SATELLITE 

- LOW ALTITUDE TETHERED SUBSATELLITE 

- SATELLITE - SATELLITE DOPPLER 






r T: fjfl 4 « 1 fjV .71 * T /»fr 4 J ), « • 


0 N S AMD UNITS 

. - cf U 

*) Xjc?Xj 

5 INDEPENDENT COMPONENTS, BY SYMMETRIES AND POISSON EQUATION 

BY MEASURING MORE THAN 5, CAN ALLOW FOR EFFECTS OF ROTATION 

UNIT OF GRAVITY : 1 GAL « 1 CM / SEC 

1 G = 10 3 GAL 

UNIT OF GRADIENT : 1 EOTVOS UNIT (EU) = 10~ 9 SEC' 2 

= 10" 9 GAL / CM 

BASELINE GRAVITY GRADIENT = 3000 EU 

$ * * 7 / VtTY * u « 2 * lO mti g U 
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IJ 
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d X, 



TETHERED SUBSATELLITE 


ADVANTAGES : 

- RESPOND TO SMALL ANOMALIES, DOWN TO SCALE ABOUT (1/ ALTITUDE) 

- STRONG SIGNAL, PROPORTIONAL TO ( 1 / ALTITUDE ) 3 

- MODERATE LENGTH MISSION 

DISADVANTAGES : 

- NOISY ACCELERATION ENVIRONMENT 

- T.S.S. NOT FULLY DEVELOPED 

















GRADIOMETRIC EFFECT 
LINEAR ACCELERATION NOISE 






REJECTED BY 

CO mm MODE ACCELERATION 


OF DYNAMIC NOISE 


TORSIONAL ACCELERATION NOISE 



PRODUCES AN ADDITIONAL GRADIENT COMPONENT: 

A = T 4 - 


ORIGINAL PAGE IS 
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SOURCES OF DYNAMIC NOISE 


- VARIABLE DRAG DUE TO ATMOSPHERIC DENSITY FLUCTUATIONS 

- TETHER OSCILLATIONS : 

- LONGITUDINAL (STRETCHING) MODES 

- LATITUDINAL MODES 

-- GROSS MOTION OF SUBSATELLITE ALONG TETHER DIRECTION ("SPRING-MASS" MODE) 

- SUBSATELLITE ROTATIONAL DYNAMICS 

- TETHER PENDULAR LIBRATIONS 

- TETHER REEL CONTROL 

- ORBITAL ALTITUDE VARIATIONS 

- MOTION ON BOARD DEPLOYED 

- MAGNETIC FIELD VARIATIONS 


it 


i 

! 
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ATMOSPHERIC DENSITY FLUCTUATIONS 


AT ALTITUDES OF INTEREST (120 - 15C KM) THE ONLY DATA IS FROM PERIGEE 
PASSAGES OF SPADES AND ATMOSPHERE EXPLORER-C SATELLITES 
~ GIVES ONLY BRIEF SAMPLES 

- VERTICAL AND HORIZONTAL CORRELATIONS CONFUSED 

- FLUCTUATIONS OF UP TO 20 % ON SCALES OF 10'S TO 100'S OF KILOMETERS 
(REBER, ET AL 1975) 

THERE IS GOOD DATA AT HIGHER ALTITUDES FROM ATMOSPHERE EXPLORER-E (300 l(M) 
AND ATMOSPHERE EXPLORER-C (250 KM) 

- LONG OBSERVATION SPANS 

- CONSTANT ALTITUDE 

- A.E.-C RESULTS (GROSS, REBER AND HUANG 1984): 

- POWER LAW SPECTRUM WITH SOME SHARP PEAKS 

- INDIVIDUAL WAVES EXTENDING ONLY A FEW PERIODS 

PRIMARY EFFECT IS HORIZONTAL DRAG ON SUBSATELLITE 

- DRAG FLUCTUATION IS PROPORTIONAL TO DENSITY FLUCTUATION 

-- BASELINE DRAG ACCELERATION AT 120 KM = 0.5 GAL = 5 x 10' 3 G 
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PROPOSED MISSIONS 


T.S.S. DEMONSTRATION FLIGHTS 
MISSION I (UPWARD): 

- STUDY SYSTEM RESPONSE USING KNOWN THRUSTER FORCES 

— COMPARE WITH THEORY, ESTIMATE DAMPING AND SYSTEM PARAMETERS 
MISSION II (DOWNWARD): 

- OBSERVE EFFECTS OF ATMOSPHERIC FLUCTUATIONS 

~ MEASURE RESPONSE AND PARAMETERS OF A DIFFERENT SYSTEM 

— POSSIBLY FLY A ONE-AXIS, NON-CRYOGENIC GRAD IOMETER TO STUDY 
G RADIOMETRIC NOISE 

DEVOTED FLIGHTS 

- FULL TENSOR GRAD IOMETER, CRYOGENIC 
— OBTAIN GEOPHYSICAL DATA 


ORIGINAL PAGE E8 
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APPENDIX D 


Detailed Plots of Simulation Results 


Figure D-l. This figure plots the acceleration perpendicular 
to the tether as a function of time, following a perturbation 
due to an atmospheric density enhancement. For details of the 
case studied, see Section 7. The two columns of numbers to the 
left of the plot are (1) the time since the beginning of the 
simulation and (2) the acceleration in gal's (cm/sec ). The plot 
begins at t = 16 seconds, since during the actual perturbation the 
(from about t=10 to t=15) the acceleration is several orders of 
magnitude greater; thus, a full span plot is simply a step function. 
Note the overall smooth appearance of the plot when compared with 
Figure D-2 below, which shows the acceleration component tangent 
to the tether. Also note the jump in acceleration near t = 420; 
an explanation for this has not been ascertained, but it is 
probably an artifact of the simulation. 
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Figure D-2. A plot similar to D-l, except that here we plot 
the acceleration component tangent to the tether. We are also 
able to plot the full span of the simulation, including the initial 
unperturbed state and during the encounter with the atmospheric 
enhancement from t = 10 to t = 16 seconds. The accelerations 
produced in this direction are several orders of magnitude 
greater than those in the orthogonal direction, and the 
maximum accelerations are comparable to the acceleration 
produced directly by the perturbing drag (0.084 gal). 
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Figure D-3, Spectral analysis of the orthogonal acceleration 
component plotted as a function of time in Figure D-l. Data 
beginning immediately after the perturbation ended, for 512 
samples, were used. Note the strong peak at low frequencies 
comprising the spring-mass mode and possibly including n 
several latitudinal tether oscillation modes (see Section 4,2.1). 
Two small peaks at somewhat higher frequencies are apparent, 
but the overall behavior at high frequencies is much smoother 
than for the tangent acceleration analyzed in Figure D-4 below. 
The logarithm of the absolute value of the spectrum is plotted. 
The columns to the left of the plot list (1) the frequency (Hz), 
(2) the associated period (seconds) and (3) the value plotted. 

As discussed in Section 6.1, the spectrum is smoothed by about 
1 frequency bin on either side of a peak because of the Hanning 
window imposed to reduce leakage of sharp spectral lines. 
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Figure D-4. A spectral analysis similar to Figure D-3 except 
that the acceleration component tangent to the tether is 
analyzed. Note the several sharp peaks out to about 0.2 Hz, 
and the jagged nature of the baseline thereafter compared with 
the baseline in D-3. The peaks correspond to the spring-mass 
mode and longitudinal tether oscillation modes as discussed in 
Section 7.1. The limited number of peaks is due to the finite 
number of masses used in the simulation. The difference in 
the baselines for the tangent and orthogonal components has 
not been analyzed. 
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